PROCRUSTES: A computer algebra package 
for post-Newtonian calculations in General 

Relativity 



Dirk Puetzfeld 

Institute of Theoretical Astrophysics, University of Oslo, P.O. Box 1029, 0315 
[ Oslo, Norway 

+-> 

O 

O 

1 Abstract 

We report on a package of routines for the computer algebra system Maple which 
supports the explicit determination of the geometric quantities, field equations, 
equations of motion, and conserved quantities of General Relativity in the post- 
Newtonian approximation. The package structure is modular and allows for an easy 
modification by the user. The set of routines can be used to verify hand calculations 
or to generate the input for further numerical investigations. 
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Program summary 



Title of the program: Procrustes 
Catalogue identifier: ADYH_vl_0 

Program obtainable from: CPC Program Library, author's webpage 



Computer for which the program is designed and others on which it has been tested: 
Computers: 

Platforms supported by the Maple computer algebra system (program was written 
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under Maple 8, but also tested with Maple 9, 9.5, 10) 

Operating systems under which the program has been tested: 
Linux, Unix, Windows XP 

Programming language used: 
Maple internal language 

Memory required to execute typical problem: 

Dependent on problem (small ~ couple of MBytes, large ~ several GBytes) 
Classification: 

1.5 Relativity and Gravitation, 5 Computer Algebra 
No. bits in a word: 

Dependent on Maple distribution (supports 32 bit and 64 bit platforms) 
No. of processors used: 1 

No. of bytes in distributed program, including test data, etc.: ~150 kbyte 
Distribution format: Compressed tar file, compressed zip archive 
Nature of the physical problem: 

The post-Newtonian approximation represents an approximative scheme frequently 
used in General Relativity in which the gravitational potential is expanded into a se- 
ries in inverse powers of the speed of light. Depending on the desired approximation 
level the field equations and equations of motion have to be determined up to given 
orders in the speed of light. This usually requires large algebraic computations due 
to the geometrical quantities entering the field equations and equations of motion. 

Method of solution: 

Automated computation using computer algebra techniques. Program has modular 
structure and only makes use of basic features of Maple to guarantee maximum 
compatibility and to allow for rapid extensions/modifications by the user. 

Typical running time: 

Dependent on problem (small ~ couple of minutes, large ~ couple of hours) 
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Restrictions on the complexity of the problem: 

Sufficient amount of memory is the limiting factor for these calculations. The struc- 
ture of the program allows one to handle large scale problems in an iterative manner 
to minimize the amount of memory required. 



1 Introduction 

Due to the structure of the field equations and the geometrical quantities en- 
tering these equations, computations in General Relativity are usually very 
tedious and prone to errors if performed by hand. This is especially true for 
calculations in which one makes use of some approximation scheme, such as se- 
ries expansions. The so-called post-Newtonian approximation, in which certain 
quantities are developed in powers of the speed of light, is an example of such 
an expansion scheme. The post-Newtonian approximation was used in many 
works to determine the equations of motion in an iterative way [18,13,11], see 
also [19] for a more comprehensive list, and plays an important role in the 
experimental verification of General Relativity [21]. 

In contrast to the efforts in the field of exact solutions, in which computer 
algebra methods are heavily used, and several specialized packages [7,2,3,1,4] 
for common computer algebra systems [5,6,9,8] exist, there do not seem to 
be many works which rely on the use of such methods when working in the 
post-Newtonian approximation. One example of a package which has been 
utilized in a post-Newtonian context can be found in [17]. Further literature 
on the use of computer algebra methods in General Relativity as well as a list 
of special purpose systems is given in [15,14]. 

In this work we present a set of routines for the computer algebra system Maple 
[5]. The routines can be used to calculate the explicit form of the field equa- 
tions, equations of motion, and conserved quantities in an automated fashion. 
The order of approximation is controlled by a set of switches and allows for 
a quick modification by the user. The output produced by the program can 
be used to support and verify hand calculations or to provide the input for 
subsequent numerical analysis. A splitting of time and resource intensive cal- 
culations into several small subproblems is supported by the modular structure 
of the routines. 

The structure of the present work is as follows. In section 2 we provide a 
general overview of the package and a listing of different types of routines. 
Subsequently, we discuss a typical application in section 3, where in we also 
give several explicit examples on the usage of different routines in the package 
and discuss their input/output. In section 4 we summarize the details of our 
package and discuss possible future extensions and modifications. Appendix 
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A contains a comprehensive directory of switches and functions available in 
Procrustes. 



2 Program structure 

2.1 General thoughts 

The set of routines which we summarize here under the name Procrustes 
emerged from the author's need to tackle a specific computation [19] within 
the post-Newtonian approximation. They were not designed to cover every 
imaginable problem in such a realm. The original task was to explicitly cal- 
culate the geometrical objects, field equations, and equations of motion for a 
given metric ansatz and to display the results in the form of a series expansion. 

Since our goal was to create a program which is easy to understand and modify, 
we did not make use of any of the existing systems for tensor manipulations. 
Consequently, there is no built-in differentiation between upper and lower 
indices of objects, just operations with lists. 

Because of the ease of use of algebraic manipulations and simplifications, we 
chose to implement our routines in the Maple [5] system. This program is 
available for many platforms and operating systems, supports large amounts 
of memory, has a large set of built in routines for symbolic computations/ 
simplifications, and supports all kinds of import and export functions. We 
only make use of very basic functions of Maple, especially its flexible list data 
structure, and intentionally use only the simplest save command, instead of 
the generation of a package repository which leads to a version dependence, 
to store our set of routines. 



2.2 Routine and program structure 

As mentioned in the introduction we are concerned with the post-Newtonian 
approximation of General Relativity. Since General Relativity represents a 
metric theory of gravitation, subsequent geometrical objects can be deduced 
from the metric. In the case of the post-Newtonian approximation the metric 
is developed into a series of inverse powers of the speed of light c starting from 

1 2 

the Newtonian limit, i.e. formally we have g a p = g5; wton + c ~ l 9ap +c~ 2 g a p 
+ . . . The strategy within the post-Newtonian approximation then typically 
consists of rewriting the field equations in a form which closely resembles the 
structures which we already know from Newtonian gravity. In addition the 
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Calculate object 


— ►> Save full object — ► 


Expand object 





Replace derivatives 



Truncate object 



Display object 



Fig. 1. Sketch of a routine which calculates and displays an object. The dashed line 
indicates that saving the object to a file is optional. 

slow-motion condition employed in the post-Newtonian approximation leads 

to a hierarchy between the different metric components when one considers a 

split into temporal and spatial components. The assumption of slow motions 

also motivates the series expansion of the velocities, which enter the energy- 

o 1 

momentum tensor, by inverse powers of the speed of light, i.e. u a =u a +c u a 

2 

+c" 2 u a + . . . We are not going into further details of the post-Newtonian 
approximation here but refer the reader to the overviews [13,16,12,20,10] and 
to the extensive literature list given in [19]. 

Our main concern, in this article, is the introduction of expansions for the 
metric g a p and velocity u a , in inverse powers of the speed of light, in all 
subsequent quantities to be calculated. In addition we must take care of the 
fact that derivatives with respect to the time coordinate increase the order of 
smallness of a term by c. 

We handle the increase of the order of smallness via derivative replacements. 
In all routines the replacement is performed after the full object has been 
calculated and only plays a role in the expression which is subsequently used 
to produce the output (see also figure 1). The display order of the output is 
directly passed to the routine calc_<obj ect name> (<display order>) . Only 
the full object remains in memory. Depending on the definition of the met- 
ric and other matter variables, which enter through energy-momentum tensor 
of the system under consideration, the user has to supply a list of substitu- 
tion rules for all occurring derivatives with respect to the time coordinate for 
these variables. This substitution list has to be provided in the global variable 
sublist. We note that for the correct operation of the substitution procedure 
it is mandatory that sublist contains a hierarchical list of time derivative 
replacements starting with the highest expected derivative. In order to help 
the user with the generation of such a substitution list the package contains a 
small routine which produces a valid list for a given set of variable names (see 
the corresponding section in the sample worksheet provided with the package). 

In addition to the routine which performs the actual calculation of an ob- 
ject, named calc_<object name>, there is also a routine, named red_<object 
name>, which performs a truncation of the object at the specified order. The 
truncation routines should be used in order to reduce the amount of memory 
needed in the calculation of subsequent objects, such as the reduction of the 
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Expand full object — > Truncate full object — >i Save reduced object 



Fig. 2. Sketch of a routine which reduces an object. The dashed line indicates that 
saving the object to a file is optional. 



Load routines 



Define switches 



Define coordinates — i 



Define metric 



Initialize substitution list 



Initialize display list 



Load / calculate / reduce / display / save objects 



Fig. 3. Preferred structure for a calculation. 

connection before the curvature is calculated. The general structure of the 
reduction routines is sketched in figure 2. The display order of the output 
is directly passed to the routine red_<object name>(<display order>). No 
derivative substitution is performed and only the truncated version of the full 
object remains in memory. Note that the reduction routines do not make use 
of the substitution list. 

The routines in the Procrustes package are simple procedures which operate 
on a set of global variables; see tables A.2-A.6 for a list of routines and table 
A. 7 for a list of variables. All routines are stored and loaded from a subdirec- 
tory defined in the global variable CCprocrustes_dir. We intentionally used 
Maple's save command to store the routines, instead of creating a package 
repository which leads to a dependence on the underlying Maple version. All 
of the routines can be found in the file procrustes_<version number> .mws. 
the user may modify this file to replace existing or add his/her own routines. 
A simple execution of this worksheet will store the routines in the package 
subdirectory. 

The recommendation for the general structure of a worksheet to be used with 
the package is sketched in figure 3. In the next section we provide some de- 
tailed comments on the input and output of the demonstration worksheet 
procrustes_demo .mws supplied together with the package. 



3 Typical application 



In this section we comment on a sample calculation with Procrustes. The 
input/output displayed here can be produced with the worksheet procrustes 
_demo.mws supplied with the package. We note that this worksheet allows one 
to quickly reproduce many of the results of the classic paper [11]. 
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3. 1 Installation 



To install the package create a directory of your choice and, depending on 
the operating system, unpack the tar or zip archive. The directory should 
contain two files procrustes_<version number>.mws (containing the source 
of the package), and pr o crust e s _demo .mws (a demonstration worksheet; ex- 
cerpts of this file are discussed in the sections below). The subdirectories 
procrustes<version number> (containing the set of routines generated from 
the source file), and testrun_procrustes<version number> (containing data 
from a testrun of the package, which is used by the worksheet to determine if 
the package is installed correctly). 

Modifications of the package can be performed by editing the source of the 
routines in procrustes_<version number> .mws. When introducing new rou- 
tines we recommend maintaining the same naming convention as introduced 
in the original file (see also the complete list of routines in tables A. 2- A. 6). 
New routines should be added to the save list at the end of the file in order to 
guarantee their storage for further usage. A simple execution of the worksheet 
then stores the modified/new routines in the package subdirectory. 

In an actual calculation we recommend to load the required set of routines 
from the package subdirectory. The worksheet procrustes_demo .mws con- 
tains, apart from some instructive sections on how to specify metrics and 
substitution lists, a sample calculation which sticks to the scheme sketched in 
figure 3. 



3.2 Coordinate and metric specification 

The coordinates to be used by Procrustes have to be specified in form of a 
list in a global variable named xx, in the first section of the worksheet this is 
done by the command: 1 

> xx : = [xO , xl , x2 , x33 ; 

Based on these coordinates the metric has to be specified in form of a list. 
The worksheet contains a section in which the contravariant as well as the 
covariant form of the metric are specified, the latter one gdd being defined as: 

> gdd:«[[l-2*U(xx[l] ,xx[2] ,xx[3] ,xx [4] )/c"2+phi (xx [1] ,xx[2] ,xx[3] ,xx [4] ) /c-4,hdd01 (xx[l] ,xx[2] ,xx[3] , xx [4] /c"3 , 

hdd02(xx[l] ,xx [2] ,xx[3] ,xx [4] /c"3 ,hdd03(xx [1] ,xx[2] ,xx[3] ,xx[4])/c"3] . 
[hdd01(xx[l] ,xx[2] ,xx[3] ,xx [4] )/c"3, (-l-2*VCxx [1] ,xx[2] ,xx[3] ,xx[4] ) /c"2) , ,0] , 
[hdd02(xx[l] ,xx[2] ,xx[3] ,xx [4] )/c"3, , (-l-2*V(xx [1] ,xxt2] ,xx[3] ,xx [4] )/c"2) ,0] , 
[hdd03(xx[l] ,xx [2] ,xx[3] ,xx [4] )/c"3,0 ,0 , (-l-2»V(xx [1] ,xx[2] ,xx[3] ,xx [4] ) /c"2)] ] ; 



In the following input by the user is marked by ">" and output by Maple is pre- 
ceeded by "<" . Three vertically aligned dots are used to indicated that we shortened 
the displayed input/output. 
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This definition corresponds to the line element 



ds 2 = g a/3 dx a dx f3 

= ( X " + ^) dx ° dx ° + ^dx°dx a - (l + ^) S ab dx a dx b , (1) 

where U, h 0a , and V represent free functions which depend on the coordi- 
nates indicated. We remind the reader that there is no built-in distinction 
between upper and lower indices, hence also the contravariant form guu of 
the metric has to be provided by hand (see sample worksheet). Coordinate as 
well as metric definitions are then saved to the files coordinates_demo .mpi , 
gdcLdemo .mpi , guu_demo.mpi for later usage in the sample calculation. 



3.3 Substitution and display list specification 



As we already have mentioned in section 2 the increase in the order of smallness 
of quantities is handled via derivative replacements. The substitution list has 
to be provided in a global variable named sublist and should contain a 
hierarchical list of replacement rules, starting with highest expected derivative 
of a variable. In the sample worksheet we have included a section in which 
a substitution list is explicitly specified. This list is subsequently used in the 
sample calculation and has the following structure: 

> sublist :=[diff (UCxO ,xl , x2 ,x3) ,x0 ,x0 ,x0 ,xO)=dif f (UCt ,xl ,x2 , x3) ,t,t,t,t)/c~4, 
diff (U(x0,xl,x2,x3),x0,x0,x0)-diff (U(t , xl ,x2 ,x3) ,t,t,t)/c"3, 
diff (U(x0,xl,x2,x3),x0,x0)-diff (UCt ,xl , x2 ,x3) ,t,t)/c"2, 
diff (U(x0,xl,x2,x3),x0)-diff (U(t,xl,x2,x3) ,t)/c, 
U(x0,xl,x2,x3)-U(t,xl,x2,x3) , 



We have added a procedure called generate_sublist to assists the user in 
the generation of a correct substitution list for all variables in a calculation. 

Apart from the substitution list, the package also makes use of a display list, 
to be defined in the global variable displist. The purpose of this list is to 
bring the output from the routines in the package into a neat form which can 
be quickly interpreted by the user. The display list does not play a role in the 
actual calculation of quantities, nor in the storage of partial results, it is only 
used when a procedure outputs an object on the screen (see figure 1). In the 
sample worksheet provided the display list is used to suppress the functional 
dependencies of the variables in the calculation, thereby reducing the length 
of the output considerably. The example for displist given therein has the 
following structure: 

> displist:=[DCt,xl,x2,x3)-D, 
V(t,xl,x2,x3)-V, 
rho Ct,xl,x2,x3)=rho, 
pp(t,xl,x2,x3)=p, 
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Depending on the needs of the user one could also use the display list to 
pre-process output which can then be directly exported into DTeX. Either 
by exclusively using the display list or by a combination of the display list 
and Maple's built-in output facility. In the next section we show some input 
and output examples from the worksheet. We caution the user that quantities 
which were reduced by the substitutions in the display list should not be used 
for any further calculations. 



3.4 Calculation and reduction of objects 



After the specification of the coordinates, metric, substitution list, and display 
list we are ready to perform a calculation with the package. In the demonstra- 
tion worksheet we provide a sample for the general structure of a worksheet for 
calculations with Procrustes. After a preamble in which we load the desired 
procedures from the package directory, and the specification of some global 
switches (see the list in table A.l), we make use of the previous definitions 
of the metric and the corresponding substitution list by loading them from 
the previously specified files and then perform some actual calculations. For 
example, the command 



> calc_uu(4) ; 



calculates and displays the contravariant four velocity up to third order in 
inverse powers of the speed of light c, 



< 222 

VU03 vuOl VU02 

+ U + + 

2 2 2 

Velocity up[, 0, ], 1 + 

2 



3 2 2 
vuOl vuOl vu03 vuOl vu02 
+ + vn01 D + 

vuOl 2 2 2 
Velocity up[, 1, ], + 

c 3 



This corresponds to the result 



u 



1 + \ (babV a V b + U 



c 2 V2 



u 



1 \5 ab v a v b + U 
- + 5 



(2) 



In addition to the actual calculation of objects it is sometimes advantageous 
to reduce the order of an object, i.e. to remove higher orders in its series 
expansion which are not needed in subsequent calculations. Such a reduction 
is supported by functions starting with red_<object name>, which exist for 
most of the objects calculable with the package, cf. table A. 5. In contrast to 
the routines which calculate objects and display them up to a given order, 
the reduction routines truncate the expanded version of the full object before 



9 



derivative substitutions are performed, see also the scheme sketched in figure 
2. As an example we calculate the connection from the sample metric, 



> calc_condndndn(4) ; 

< d 

— U 
dt 

Connection (condndndn) [ , 0, 0, 0, ], - 

3 



Note that the order passed to the calculation command only influences up to 
the order objects are displayed, not the actual form of the object which is kept 
in memory for further calculations. We can display the full object by directly 
accessing the corresponding list element, which in this case is: 



> condndndn [1 , 1 , 1] ; 



d d 

— DCxO, xl, x2, x3) phi(x0, xl , x2, x3) 

dxO dxO 

+ 1/2 

2 4 



Hence, we see that the full object is still in memory. Only the reduce command 
truncates the full object and thereby affects subsequent calculations, as we can 
see from the following command sequence: 



> red_condndndn(4) ; condndndn [1 ,1,1]; 



d 

DCxO, xl, x2, x3) 

dxO 



2 



This also shows that no derivative replacements are performed during the 
reduction of full objects. The maximum order up to which objects get trun- 
cated is controlled via the global variable CCmaxorder (cf. table A.l). The 
user should make sure that the order specified in this variable is high enough 
to guarantee that all unwanted higher order terms are removed. If the user 
wants to display a single component of an object including the derivative re- 
placements the command disp_component, cf. table A. 2, should be used. For 
the connection component from the previous examples this command yields 



> disp_component C ' condndndn [1 , 1 , 1] ' ,4) ; 



d 

— U 
dt 

condndndn [1, 1, 1], 

3 



Finally, we note that the reduction commands should be used with great care, 
since they directly influence the calculation of subsequent objects if these 
depend on the object being reduced. 
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3.5 Saving and loading components 



Time and memory intensive calculations are supported by the package via an 
integrated save function, which is controlled by the global switch CCoutputf s 
(see table A.l). If this switch is set to the value 2, all functions produce a 
file with reusable output. The calculation routines produce output files with 
names <ob j ect/component name> . full . mpi. Files with the names <ob j ect/ 
component name> . red. mpi are generated by the reduction routines. Output 
files are generated for every component of an indexed object and there are 
routines for loading the entire object as well as only single components. The 
sample worksheet contains a section in which the previously computed results 
are loaded from files. For all objects there exists a corresponding load routine 
load_<object name> (full/reduced switch) which reads the complete ob- 
ject from files, the switch passed to the load command determines whether files 
with the suffix .full. mpi or .red. mpi are used (see table A. 6 for a complete 
list). 

Components of objects can be loaded from the existing files by using the 
load_component command. In contrast to the load_<object name> command 
the load_component command returns the specified component of an object 
but does not automatically assign it to the corresponding variable name. Here 
are some examples for objects with one and two indices and a scalar quantity: 

> load_ component ( 'uu' , [1] , 1) ; load_ component ( 'tupup ' , [1 ,2] , 2) ; load_component ( ' thetascalar ' , [] ,2) ; 

< "uu_l . full .mpi" 
"tupup_l_2.red.mpi" 
"thetascalar. red. mpi" 

As noted before the above commands do not automatically perform assign- 
ments to the corresponding variable names. This can be easily done by hand 
as follows 

> for aa from 1 to 4 do test_uu [aa] : =load_ component ( 'uu' , [aa] , 2) od: 

< "uu_l .red. mpi" 
"uu_2.red.mpi" 
"uu_3.red.mpi" 
"uu_4.red.mpi" 

or by using the predefined load command, which loads all components of uu 
from the corresponding files and assigns them. We can verify the equivalence 
of the last command sequence and the load command by 

> load_uu(2) ; 

< "uu_l .red. mpi" 
"uu_2.red.mpi" 
"uu_3.red.mpi" 
"uu_4.red.mpi" 

> seq(is(test_uu[aa]=uu[aa] ) ,aa=l . .4) ; 

< true , true , true , true 

More examples can be found in the sample worksheet. The complete list of 
predefined load commands is given in table A. 6. 
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3.6 Testrun 



After the installation it is advisable to perform a quick testrun with the pack- 
age. This can easily be done by executing the worksheet supplied with the 
package. After performing the sample calculations (which should not take 
longer than a couple of minutes on a modern desktop machine 2 ) the work- 
sheet automatically compares its results with some test data files which are 
supplied together with the package (as noted before the data for comparison is 
located in the directory testrun_procrustes<version number>). All of the 
current features of the package are demonstrated in the sample worksheet. We 
note that the calculation in the sample worksheet reproduces many results of 
the classic work [11] and is therefore a good example for the effectiveness of 
the package. 



4 Conclusions and outlook 



We have introduced the main features of the Procrustes package and explained 
how it can be used to tackle computations within the post-Newtonian approx- 
imation of General Relativity. The package is versatile and should allow the 
user to easily make modifications. Possible future developments of the package 
strongly depend on the scientific problems in which the author will be involved. 
Planned future developments include: (i) An automated order control for all 
quantities entering a specific post-Newtonian calculation. This feature would 
further reduce the memory requirements of the program and would alleviate 
the user from making choices for the correct post-Newtonian orders of single 
quantities, (ii) Automated simplification schemes for the output. The goal is 
to introduce some extra routines which match certain patterns in the output 
and thereby reduce the need for manual calculations even further, (iii) Ex- 
tension of the package to incorporate non-Riemannian theories of gravitation. 
The complexity of the geometric quantities in such theories renders virtually 
any hand calculation unfeasible. From a physical point of view this is one 
of the most pressing steps, since it would allow us to establish a systematic 
post-Newtonian framework for alternative gravitational theories. 



2 With Maple 8 it takes about 35 seconds on a Pentium M (1.7 GHz) system and 
approximately 320 MBytes of main memory. 
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Table A.l 

List of global switches. 



Name 


Description 


CCmaxorder 


Order up to which expanded objects get truncated. The 




standard setting is 200, meaning that the reduction of ob- 




jects is performed starting from the order passed to the rou- 




tine up to c -200 


CCoutputf s 


Switch which toggles the output to files. The standard set- 




ting is 2, meaning that all routines will produce files with 




their output 


CCprocrustes_dir 


Directory in which the routines of the package are stored, 




the standard setting is " . /procrustes<versionumber>/" 


CCtaylor 


Order variable used in Taylor approximations, the standard 




setting is 3 
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A Directory of functions and names 

In this section we provide the names and a description of predefined switches 
A.l, a list of routines and their dependencies A.2-A.6, and a list of global 
variables A. 7 in table form. 
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Table A.3 

List of predefined functions which calculate an object and display it up to the given 
<display order>. 



Name 


Dependencies 


Description 


calc_uu(<di splay order>) 


gdd, fapproxQ 


Contravariant four velocity u a = dx a /ds 


calc_ud(< display order>) 


gdd , uu 


Covariant four velocity u a = g a p u ^ 


calc_tupup(<di splay order >) 


uu , guu 


Contravariant perfect fluid energy- 
momentum tensor T a @ 


calc_tdndn(<di splay order >) 


ud , gdd 


Covariant perfect fluid energy-momentum 
tensor T a p = (pc 2 + lip + p) u a Uf3 — pg a p 


calc_tupdn(<di splay order >) 


gdd, tupup 


Mixed perfect fluid energy-momentum 

IjClloUl -L Q 


calc_rhs (<display order>) 


gdd, tdndn, tupdn 


Right-hand side of the field equations 
RHS a/ g = (T a p — ^T 7 7 (/ Q/ g) 


calc_condndndn(<di splay order >) 


gdd 


Connection T^ af3 = \ (g ja ,p + gp 7 , a 

— do 1 


calc_conupdndn(<di splay order >) 


guu , condndndn 


Connection 1 a/3 — g ' 1 S \ a/3 


calc_ctupup(<di splay order>) 


tupup , conupdndn 


Covariant derivative of the EM tensor 
1 ;i 






Equations of motion T a ^ -p — 


calc_ricdndn(<display order>) 


guu , gdd , condndndn 


Covariant Ricci tensor R^ v = 

_!_ // ^ ( ri n _i_ // i< n / j i ^ _|_ 

2i/ \yfJ,v,af3 T gaf3,fiv yau,fj,p y^P,au)^ 

5 a ^5 75 (r 5 | j[1I/ r 7 | a/ g -T^pT^v) 


calc_ricupdn(<display order>) 


guu, ricdndn 


Mixed Ricci tensor R a p 


calc_ricscalar (<di splay order >) 


ricupdn 


Mixed Ricci scalar R = R a a 


calc_riemupdndndri(<di splay order>) 


conupdndn 


Riemann curvature tensor R a p^v = 

-pa i -pa -pi pa _i_ -pi -pa 
1 P^,v 1 Pv,fjL 1 Pn iv "•" 1 Pu 1 7M 


calc_riemdndndndnC<di splay order>) 


gdd, riemupdndn 


Riemann curvature tensor R a p^,v = 

<7a7 R? p^v 


calc_weyldndndndri(<di splay order>) 


gdd, ricdndn, 
ricscalar , 
r i emdndndndn 


Weyl curvature tensor C a p^ u = Rap^u — 
\{gav Bp^ + gp^ B au — g a/x Bp v —gp v B a/X ) — 

^R(g a ugpfi - ga^gpv) with B a p = R a p - 
jgapR 


calc_detguu(<di splay order>) 


guu 


Metric determinant det(g a ^) 


calc_detgdd(<di splay order>) 


gdd 


Metric determinant det(g a p) 
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Table A.4 

List of predefined functions which calculate an object and display it up to the given 
<display order>. 



Name 



Dependencies Description 



calc_pemtupup(<di splay order >) 



guu , conupdndn 



calc_thetaupup(<di splay order >) 



detgdd, tupup, 
pemtupup 



calc_thetaupuppdn(<di splay order>) thetaupup 



calc_harm_gc (<display order>) 



guu, detgdd, 
gapproxQ 



calc_harm_gc_alt (<display order>) guu, conupdndn 



calc_pn_gc (<display order>) 

calc_cuu(<di splay order>) 

cal c_ cud (<di splay order>) 

calc_thetascalar (<display order>) 
cal c_acd(<di splay order>) 
cal c_omegadd(< display order>) 
calc_thetadd(<display order>) 

calc_sigmadd(<display order>) 
check_velnorm(<di splay order >) 



Landau-Lifshitz 
pseudotensor 



energy- moment um 

„ua „8S (tv rn i rf pk -pu pb 

9 9 {^aK L /3S + L /3S L a K ~ 1 Sk L a/3 ~ 

L al3 L 8K> +99 [LccKl-ps i L /3S L aK 

HM* - rj,a) + 9 aP 9 SK K s ^ K - 

Landau-Lifshitz energy-momentum com- 
plex 6 Q/3 = - det(g al3 )[T a P + t af3 } 



"p5 ~pK 

1 as 1 /36 



Partial derivative of the LL EM complex 



Harmonic gauge condition (y/—gg' 



af3\ 



gdd 



uu, conupdndn 



ud, conupdndn 



uu , cud 



ud , cud , acd 



ud , cud , acd 



gdd , ud , thetadd , 
thetascalar 



uu, ud 



Harmonic gauge condition <7 a ^r^g (alter- 
native form) 

Standard post-Newtonian gauge condition 
(S ab (9oa, b -k9ab, ) = 0, 5 hc g ab , c -\{5 bc g bc - 

900),a = 0) 

Covariant derivative of the four velocity 
Covariant derivative of the four velocity 

U a -/3 

Expansion scalar 9 



u 



che ck_metric id (<di splay order >) guu, gdd 

che ck_de compos it ion (<di splay order >) 



omegadd , s igmadd , 
thetascalar, gdd, 
ud , acd , cud 



Acceleration a a = u a -pu 



Rotation u a/3 = U[ a .p\ + U[ a a^ 

Symmetric part of the derivative projec- 
tion 6 af3 = U( a;y9 ) - u (a aa) 

Shear a a/3 = 6 a p - \9(g a p - u a up) 

Checks fulfillment of the velocity normal- 
ization 

Checks fulfillment of the metric identity 

Checks whether the decomposition n Q;/ g = 
u a (i + cr a /3 + \{9af3 - u a U/3)0 + a a U/3 is con- 
sistent 
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Table A.5 

List of predefined functions which truncate an object at the given <truncation 
order>. 



Name 


Dependencies Object 


red_uu(<truncation order>) 


uu 


Four velocity u a 


red_ud(<truncation order>) 


ud 


Four velocity u a 


red_tupup(<truncation order>) 


tupup 


EM tensor T a/3 


red_tdndn(<truncation order>) 


tdndn 


EM tensor T a/3 


red_tupdn(<truncation order>) 


tupdn 


EM tensor T a p 


red_rhs (<truncation order>) 


rhsf eqsdndn 


Right-hand side of FEQs RHS a/3 


red_condndndn(<truncation order>) 


condndndn 


Connection 


red_conupdndn(<truncation order>) 


conupdndn 


Connection 


red_ctupup(<truncation order >) 


c tupup 


Covariant deriv. of the EM tensor T al3 - n 


red_eom(<truncation order >) 


eom 


EOMs T^.p = 


red_ricdndn(<truncation order >) 


ricdndn 


Covariant Ricci tensor R a p 


red_ricupdn(<truncation order >) 


ricupdn 


Mixed Ricci tensor R a p 


red_ricscalar (<truncation order>) 


ricscalar 


Ricci scalar R a a 


red_riemupdndndn(<truncation order >) 


riemupdndn 


Riemann curvature tensor R a 


red_riemdndndndn(<truncation order >) 


riemdndndn 


Riemann curvature tensor Rap^u 


red_weyldndndndn(<truncation order >) 


weyldndndndn 


Weyl curvature tensor C a p^ y 


red_detguu(<truncation order >) 


detguu 


Metric determinant det(g a @) 


red_detgdd(<truncation order >) 


detgdd 


Metric determinant det{g a p) 


red_pemtupup(<truncation order >) 


pemtupup 


LL EMPT t a ? 


red_thetaupup(<truncation order>) 


thetaupup 


LL EM complex 9 Q/3 


red_thetaupuppdn(<truncation order >) 


thetaupuppdn 


Partial deriv. of the LL EM complex ,/3 


red_cuu(<truncation order >) 


cuu 


Covariant derivative u a -p 


red_cud(<truncation order >) 


cud 


Covariant derivative u a -p 


red_thetascalar (<truncation order >) 


thetascalar 


Expansion scalar 


red_acd(<truncation order >) 


acd 


Acceleration a a 


red_omegadd(<truncation order >) 


omegadd 


Rotation u a p 


red_thetadd(<truncation order >) 


the t add 


Symm. part of the deriv. projection 6 a p 


red_sigmadd(<truncation order >) 


sigmadd 


Shear a a p 
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Table A.6 

List of predefined functions which load objects from files. The parameter 
<full/reduced expression> passed to the function determines whether the 
full or reduced object is loaded (1 = <object name> . f ull .mpi , 2 = <object 
name>.red.mpi). All functions depend on the load_component () routine. 



Name 


Object 


loacLuu(<full/r educed expression>) 


Four velocity u a 


load_ud(<f ull/r educed expression>) 


Four velocity u a 


load_tupup(<f ull/reduced expression>) 


EM tensor T a/3 


load_tdndn(<f ull/r educed expression>) 


EM tensor T af3 


load_tupdn(<f ull/reduced expression>) 


EM tensor T a p 


load_rhs (<f ull/reduced expression>) 


Right-hand side of FEQs RHS Q/3 


load_condndndn(<f ull/r educed express ion>) 


Connection T^ a p 


load_conupdndn(<f ull/reduced express ion>) 


Connection Via 


load_ctupup(<full /reduced expression>) 


Covariant derivative of the EM tensor T a/3 ;7 


load_eom(<f ull/reduced expression>) 


EOMs T a P. p = 


loader icdndn(<f ull/reduced expression>) 


Covariant Ricci tensor R a p 


load_ricupdn(<f ull/reduced expression>) 


Mixed Ricci tensor R a p 


load_ricscalar (<f ull/reduced expression>) 


Ricci scalar R a p 


load_riemupdndndn(<f ull/reduced expression>) 


Ricmann curvature tensor R a /3 fll/ 


load_riemdndndndn(<f ull/reduced expression>) 


Riemann curvature tensor R a /3nv 


load_weyldndndndn(<f ull/reduced expression>) 


Weyl curvature tensor Cap^v 


load_detguu(<f ull/reduced expression>) 


Metric determinant det(g a/3 ) 


load_detgdd(<f ull/reduced expression>) 


Metric determinant det{g a p) 


load_pemtupup(<f ull/reduced expression>) 


LL EMPT t a P 


load_thetaupup(<f ull/reduced expression>) 


LL EM complex 6 a/3 


load_thetaupuppdn(<f ull/reduced expression>) 


Partial derivative of the LL EM complex <d a @ p 


load_cuu(<f ull/reduced expression>) 


Covariant velocity derivative u a -p 


load_cud(<f ull/reduced expression>) 


Covariant velocity derivative u a -p 


load_thetascalar(<f ull/reduced expression>) 


Expansion scalar 9 


load_acd(<f ull/reduced expression>) 


Acceleration a a 


load_omegadd(<f ull/reduced expression>) 


Rotation u a p 


load_thetadd(<f ull/reduced expression>) 


Symmetric part of the derivative projection 9 a p 


load_sigmadd(<f ull/reduced expression>) 


Shear a a p 
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Table A.7 

List of global variables 



Name 



Description 



Kinematic acceleration a a 
Connection T^ a p, 

Covariant derivative of the EM tensor T a/3 -j 
Covariant derivative of the four velocity u a - t p, u a -p 
Metric determinant det(g a p), det(g a/3 ) 
Substitution list for replacements in output 
Equations of motion T a/3 .p 
Covariant g a p / contravariant g alS metric 
Harmonic gauge condition 
Kinematic rotation u> a p 
Landau-Lifshitz EM pseudotensor t a/3 
Standard post-Newtonian gauge condition 
Ricci tensor R a p, R a p 
Ricci scalar R 

Riemann curvature tensor R a p^s, R a p^8 

Right hand side of the field equations — ^ (T a p — \T^ 1 g a p ) j 



acd 

condndndn , conupdndn 

ctupup 

cuu , cud 

detgdd, detguu 

displist 

eom 

gdd, guu 

harm_gc , harm_gc_alt 

omegadd 

pemtupup 

pn_gc 

ricdndn , ricupdn 
ricscalar 

r i emupdndndn , r i emdndndndn 

rhsf eqsdndn 

routine_list 

sigmadd 
sublist 

tdndn , tupup , tupdn 
thetadd 

thetaupup, thetaupuppdn 



ud, uu 



weyldndndndn 



List containing the names of routines to be loaded from rou- 
tine directory 

Kinematic shear a a p 

Substitution list for derivative replacements 

Energy- momentum tensor T a p, T a/3 , T a p 

Kinematic quantity 9 a p, symmetric part of the velocity 
derivative projection 

Landau-Lifshitz EM complex Q/3 and its partial derivative 
Kinematic expansion 9 

Covariant u a / contravariant u a four velocity 
Velocity product u a u a 
Weyl curvature tensor C a p^s 
List with coordinates 



19 



[16] L. Infeld and J. Plebanski. Motion and Relativity. Pergamon Press, New York, 
1960. 

[17] S. A. Klioner. New system for indical computation and its applications in 
gravitational physics. Comp. Phys. Comm., 115:231, 1998. 

[18] H.A. Lorentz and J. Droste. The motion of a system of bodies under the 
influence of their mutual attraction, according to Einstein's theory. Versl. K. 
Akad. Wet. Amsterdam, 26:392, 1917. 

[19] D. Puetzfeld. The cosmological post-Newtonian equations of hydrodynamics in 
General Relativity. Submitted, 2006. 

[20] M.H. Soffel. Relativity in astronomy, celestial mechanics, and geodesy. Springer, 
Berlin, 1989. 

[21] CM. Will. Theory and experiment in gravitational physics. Cambridge 
University Press, Cambridge, 1993. 



20 



